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ABSTRACT 


Heat-transfer  coefficients  have  been  determined 
for  a  crossflow  tube  bank  by  the  use  of  heat-and  mass- 

* 

"'w  '/transfer  techniques  for  the  Reynolds -number  range  of 

35,000  to  80,000.  The  average  heat-transfer  coefficients 
agree  favorably  with  extrapolations  of  existing  data.  The 
local  coefficient  variations  from  the  average  were  as  great 
as  +80  per  cent  and  -55  per  cent. 
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AVERAGE  AND  LOCAL  HEAT  TRANSFER  FOR  CROSSFLOW  THROUGH  A  TUBE  BANK 
By  R.A.  DeBortoli,  R.E.  Grimble,  and  J.E.  Zerbe 

NOMENCLATURE 

The  following  nomenclature  Is  used  In  the  paper: 

Roman  Letter  Symbols 

Cc  «*  specific  heat  capacity  of  test-rod  material 

Cp  ®  specific-heat  capacity  of  coolant 

D  *»  diameter  of  test  rod 
h  =  convection  heat -transfer  coefficient 
hj.  -  radiation  heat-transfer  coefficient 

h  •  total  heat-transfer  coefficient 
t 

K  «  molar  mass -transfer  coefficient  of  naphthalene 
m  thermal  conductivity  of  coolant 

k  »  molar  mass -diffusion  coefficient  of  naphthalene 
d 

L  **  effective  length  of  test  rod 

M  ■  molecular  weight  of  air 
a 

M  «  molecular  weight  of  naphthalene 
n 

P  «  total  static  pressure 

p  ■»  vapor  pressure  of  naphthalene 

q*  «  flux  of  radiant  energy 
r 

T  «  absolute  temperature  of  rod 
T&  »  absolute  temperature  of  air 

t  «  time 


Accesion  For 

NTIS 

DTIC 

Unanrn 

Justific 

CRA&I  V 

TAB  m 

junced  ™ 

ation 

By . 

Distribution/ 

Availability  Codes 

Dist 

H 

a* 

Avail 

Spc 

ind  /  or 
cial 

v 


maximum  velocity  of  air 


-  2 


W  «  mass  of  test  rod 
o 


y  =  depth  of  surface  removal  from  naphthalene  rod 

Greek  Letter  Symbols 

€:  =  emissivity  of  rod  surface  for  radiation 
9  =  total  time  for  testing 

JJ~  =  dynamic  viscosity  of  air 
-  *  density  of  air 


l°n 


density  of  naphthalene 


(f  =  Stefan  Boltzmann  constant  for  radiation 


Dimensionless  Parameters 


heat  transfer  factor 


CP  Pa  V 


KP 

3  =  mass  transfer  factor  ■-  -  — 

■  Pa  y/»a  l  Pa  *d 


hD 

Nu  =  Nusselt  number  - — 

k 

a 


Pr  m  prandtl  number 


2/3 


2/3 


Re  •>  Reynolds  number 


Dv 


INTRODUCTION 

Heat-transfer  data  for  crossflow  through  tube  banks  are  abundantly  reported  In 
the  literature  for  the  lower  Reynolds  numbers.  A  range  of  Prandtl  numbers  and  a  variety 
of  geometric  arrays  have  been  covered,  and  correlations  have  been  reported.  There  is, 
for  example,  the  Colburn  equation  {1) ,\J/ 

Nu  -  0.33  (Pr)1^  (Re) 


^Underlined  numbers  in  parentheses  refer  to  the  Bibliography  at  the  end  of  the  paper. 
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that  is  based  on  data  taken  in  the  interval  from  Re  =  2000  to  Re  =  32,000.  Available 
data  become  more  sparse  as  higher  Reynolds  numbers  are  considered .  The  paper  by 
Sheehan,  et  al.  (_2),  who  used  water  as  a  coolant  at  a  temperature  of  360  F,  reports 
the  only  data  known  by  the  authors  to  have  been  taken  in  the  range  of  Reynolds  numbers 
from  105  to  10^. 

It  was  felt  that  additional  data  for  the  high  and  intermediate  Reynolds  numbers 
would  be  of  value.  The  range  covered  in  the  work  reported  here  is  between  Re  «  35,000 
and  Re  *  80,000.  Air  was  used  as  the  flow  medium  and  the  prandtl  was  taken  as  O.697. 

The  geometry  was  that  of  a  bank  of  staggered  rods  so  spaced  that  the  centers  of  any 
three  mutually  adjacent  rods  lay  at  the  vertices  of  an  equilateral  triangle  the  length 
of  whose  sides  was  1.30  times  the  rod  diameter. 


AVERAGE  HEAT -TRANSFER  COEFFICIENTS 

Baals  of  the  Method  of  Test 

The  transient  technique,  reported  by  Kays,  et  al.  (j5)  was  used  in  the  average 
heat-transfer  measurements  reported  here.  With  this  method,  if  one  rod  in  a  bank  of  rods 
is  preheated  and  is  then  cooled  by  air  flow,  its  temperature  changes  according  to  the 
equation 

d<T  -  Ta>  -  -  V51  (T  -  T„ ) 

d  t  1  q 

"c^c 

Here  it  is  assumed  that  h  and  T„  are  constant  in  time  and  that  the  rod  is  uniform  in  tem- 

t  a 

perature.  For  convenience  in  computation,  this  equation  may  be  written 


-2-  ln(T«  -  T£)  -  -  W  DL  . (1) 

d  t  w  c 

c  c 


where  T’  -  T’  is  any  quantity  proportional  to  temperature  difference  T  -  Ta  (for  example, 
a  thermocouple  potential).  Thus,  a  plot  of  ln(T'  -  T£)  against  t  will  result  in  a 
straight  line  of  slope 


V 

V 


h^  it  DL 


Apparatus  and  Test  Procedure 

The  basic  model  for  the  tests  is  shown  in  Fig.  1.  It  consisted  of  nine  rows  of 
plexiglass  rods,  1  in.  dlam.  The  rods  were  staggered  so  as  to  form  an  equilateral  tri¬ 
angular  array.  The  transverse  pitch  was  1.30  in.  and  the  longitudinal  pitch  was  1.125  in. 
The  rods  were  held  in  a  rectangular  plexiglass  channel  section.  The  inside  cross-section* 

al  dimensions  were  5  5/8  in.  on  the  side  parallel  to  the  rod  axis  and  6  3/8  in.  on  the 
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side  normal  to  the  rod  axis.  Each  odd-numbered  row  contained  five  rods,  and  each 
of  the  even-numbered  rows  had  four  whole  rods  and  a  partial  rod  at  either  end.  Any 
one  of  the  three  rods  at  the  centers  of  rows  3,  5,  and  7,  counting  In  the  direction  of 
flow,  could  be  removed  through  the  sides  of  the  channel  section  and  replaced  by  an 
Instrumented,  preheated  copper  rod  that  served  as  a  thermal  capacitor. 

Two  different  copper  test  rods  were  used.  One  consisted  of  a  0.510-lb  cy¬ 
linder,  4  5/8  In.  In  length,  1.00  In.  0D,  and  3/4  In.  ID.  The  second  capacitor  weighed 
0.864  lb  and  had  the  same  dimensions  as  the  first,  except  that  the  Inside  diameter  was 
1/2  In.  A  drawing  of  this  rod  Is  shown  in  Pig.  2.  in  order  to  eliminate  wall  effects, 
the  copper  tests  rods  were  1  in.  shorter  than  the  dummy  rods  so  as  not  to  span  the  entire 
distance  between  the  walls  of  the  test  section.  Extra  length  for  support  and  preser¬ 
vation  of  flow  conditions  was  effected  by  mounting  the  rods  between  Teflon  end  pieces  of 
the  same  outside  diameter.  Thermal  contact  between  the  copper  and  the  Teflon  was  kept 
to  a  minimum  by  recessing  grooves  in  the  Teflon. 

A  rapid-response  thermocouple  Junction  was  made  by  drilling  two  small  holes, 
close  together,  through  the  wall  of  the  cylinder  and  soldering  an  Iron  wire  In  one  hole 
and  a  const ant an  wire  in  the  other.  The  cold  junction  was  placed  in  the  air  stream  so 
that  the  output  potential  of  the  thermocouple  indicated  the  difference  in  temperature 
between  the  rod  and  the  air.  The  thermocouple  leads  were  connected  to  each  of  the  six¬ 
teen  pairs  of  terminals  of  a  Leeds  and  Northrup  Speedomax  l6-point  recorder.  The  rod 
temperature  was  thus  recorded  at  3-sec  intervals.  The  recorder  was  set  on  a  range  of  0 
to  2.5  mv,  corresponding  to  a  temperature  range  of  approximately  80  deg  F. 

Air  flow  was  maintained  by  a  Buffalo  Porge  Blower  that  was  driven  by  a  variable- 
speed  motor.  During  runs  at  the  higher  flow  rates,  the  air  passed  through  a  rectangular 
exhaust  section  and  flow  rate  was  determined  by  traverses  with  a  Hastings  precision  air 
meter.  The  flow  rate  was  determined  at  the  lower  flow  rates  by  forcing  the  air  through 
a  2.5-ln.  nozzle  and  measuring  velocity  by  a  pitot-static  primary  element  connected  to  a 
micromanometer . 

To  run  a  test,  the  blower  speed  was  adjusted  to  the  desired  flow  rate  while  the 
copper  test  rod  was  heated  on  a  hot  plate  to  about  250  F.  The  test  rod  was  then  sub¬ 
stituted  for  one  of  the  dummy  rods.  The  recorder  began  charting  the  temperature -time 
history  of  the  rod  as  soon  as  the  thermocouple  output  fell  within  the  recorded  range.  As 
the  test  rod  cooled,  the  flow-rate  data  were  taken,  the  air-stream  temperature  was  measured 
with  a  mercury  thermometer,  the  pressure  In  the  test  section  above  atmospheric  was  measured 
by  means  of  a  water  manometer,  and  the  barometer  reading  was  recorded. 

Evaluation  of  Data 


The  excess  of  rod  temperature  over  air-stream  temperature  was  read,  point  by  point, 
and  replotted  as  a  function  of  time  on  two-cycle,  semi-log  co-ordinate  paper  and  the  most 
representative  straight  line  was  drawn,  as  shown  In  Fig.  3.  The  slope  was  then  calculated 
from  two  widely  separated  points  read  from  the  straight  line.  The  total  heat-transfer 
coefficient  was  then  calculated  from  the  equation 


In  (T2  -  Ta)  -  In  (Tx  -  Ta) 


t  -  t 
2  1 


ht  it  DL 


Vo 


(2) 
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Correction  for  Radiation 


The  heat-transfer  coefficient  as  calculated  by  the  foregoing  procedure  Includes 
a  small  radiation  effect.  The  net  flux  of  radiant  energy  from  the  surface  of  the  test 
rod  Is  given  by  the  equation 


qj  -  86  (T4  -  Ta4)  -8S(T  -  Ta)  (T  +  Ta)  (T2  +  Ta2) 


assuming  that  the  surrounding  dummy  rods  comprise  a  black  body  at  air-stream  temperature. 
Then 


hp  -  -qf._ . -  =  <$6(T  +  Ta)(T2  +  Ta2)  ^  4  8^Ta  2 . (3) 

T  -  Ta 

where  TaVg  la  the  average  absolute  temperature  of  the  rod  during  a  run. 

The  surface  of  the  test  rod  changed  by  oxidation  as  it  was  exposed  to  the  air  and 
by  cleaning  between  tests,  but  €  -  0.60  was  taken  as  a  reasonable  average  value  for  cal¬ 
culating  hj, .  The  measured  heat-transfer  coefficient  was  then  corrected  by  the  relation 
h  •  ht  -  hp.  The  size  of  this  oorreetlon  was  from  1.5  per  cent  to  3.5  per  cent  of  the 
total  value. 

Por  generality,  the  heat-transfer  coefficients  were  plotted,  as  Nu/Pr1^  against 
flow  converted  to  Re,  on  logarith  co-crdlnate  paper. 

Estimate  of  Error 

An  attempt  has  been  made  to  estimate  the  magnitude  of  errors  that  might  appear  In 
the  final  results  due  to  unsuspected  causes.  Plrst  we  observe  that  the  final  heat-transfer 
coefficient  Is  given  lr  terms  of  measured  quantltltes  by  the  formula 


h 


wcc0  m  (Tg  -  Ta)  -  m  (Tl  -  Ta)  _  ^ 

IT  DL  t2  -  tj 


.(4) 


It  was  conservatively  estimated  that  the  odds  were  twenty  to  one  that  the  error  In  these 
measured  quantities  fell  within  the  following  intervals: 

W0  +0.01  lb 

C0  +0.005  Btu/lb  deg  P 

D  +0.006  In. 

L  +0.010  In. 


J 
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T-,  -  Ta  +0.02  mv 

T2  -  Ta  +0.02  mv 

tg  -  t^  +5.0  sec 

hr  +0 .4  Btu/hr  sq  ft  deg  P 

Following  the  procedure  described  by  Kline  and  McClintock  (4_) ,  it  was  calculated  for  a 
typical  run  that  the  odds  were  twenty  to  one  (95  par  cent  confidence  limit)  that  the 
fractional  error  in  the  heat-transfer  coefficient  fell  within  the  interval  +7>5  per  cent. 

Discussion  of  Results 

The  published  work  of  other  experiments  shows  that  the  heat-transfer  coefficient  in  a 
staggered  tube  bank  Increases  from  the  first  row  to  the  fourth,  reaching  essentially  its 
asymptotic  value  at  the  fourth  row  . 

For  this  reason,  the  results  obtained  with  the  test  rod  in  Row  3  were  plotted 
separately,  as  Nu/Pr1/^  versus  Re,  in  Fig.  4.  In  the  lower  Reynolds  nmriber  range  these  results 
are  compared  with  the  Colburn  Equation  (1)  into  which  has  been  Inserted  the  factor  O.93,  recom¬ 
mended  by  Kays  and  Lo  (£) .  In  the  higher  range  of  Reynolds  numbers,  the  results  are  compared  with 
the  correlation  of  Sheehan,  et  al.  ( 2 )  who  found  the  asymptotic  coefficient  to  have  been  reached 
by  the  third  row.  The  Sheehan  correlation  is  given  as 

h  m  0.185  (Re)0'®  Btu/hr  sq  ft  deg  F 

This  equation  becomes 

Nu  -  0.0331  (Re)0'8  (Pr)1/3 


assuming  that  Nusselt  number  varies  with  prandtl  number  to  the  one-third  power. 

Heat-transfer  coefficients  measured  in  the  fifth  and  seventh  rows  are  plotted  in 
Fig.  5  and  are  correlated  by  the  equation 


Nu  =  0.126  (Re)0'692  (Pr)1,/3 


(5) 


Again,  the  results  are  compared  with  the  correlation  of  Sheehan,  et  al.  and  with  the  Colburn 
equation.  These  values  are  seen  to  be  slightly  higher  than  would  be  predicted  by  either 
the  Colburn  equation  or  the  Sheehan  equation.  This  tendency  appears  reasonable,  however, 
since  a  smooth  Interpolation  between  the  two  curves  must  fall  above  the  extrapolation  of 
either . 


The  third-row  results  are  not  significantly  different  from  the  fifth-and  seventh- 
row  results.  This  is  in  agreement  with  the  findings  of  Sheehan,  et  al.,  but  not  with  the 
data  of  Kays  and  Lo. 
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VARIATION  OF  LOCAL  HEAT-TRANSFER  COEFFICIENTS 


Basis  of  the  Test 


Variations  in  heat  transfer  around  the  circumference  of  a  rod  in  crossflow 
were  determined  from  mass-transfer  data  using  the  mass  transfer  ■  heat  transfer  analogy. 
The  mass  transfer  J -factor  is  defined  to  be 


KP 


Pa  V/Ma 


.(6) 


and  is  shown  by  Chilton  and  Colburn  (5)  to  be  equal  under  the  proper  conditions  to  the 
heat-transfer  j -factor. 


J 


h 


h 

cp  Pav 


V±) 


2/3 


k 

a 


(7) 


The  conditions  under  which  this  equality  holds  are  closely  approximated  in  the  sublimation 
of  naphthalene  into  an  air  stream  at  room  temperature . 

Description  of  Apparatus  and  Test  procedure 

The  same  basic  rod  assembly  and  the  same  flow  equipment  were  used  in  these  tests 
as  were  used  for  the  transient  measurements  of  average  heat-transfer  coeflcients,  with  the 
exception  that  a  naphthalene  rod  replaced  the  heated  copper  rod. 

The  test  rods  were  cast  from  sublimed  flake  naphthalene  in  a  split  brass  mold  of 
inside  diameter  1.004  in.  After  removal  from  the  mold,  the  rod  was  cut  to  the  correct 
length  and  the  regions  next  to  the  ends  were  wrapped  with  Scotch  celulose  tape. 

When  the  desired  flow  rate  was  obtained  by  adjusting  the  blower  speed,  a  plexi¬ 
glass  rod  was  removed  and  replaced  by  a  naphthalene  rod.  The  rod  was  oriented  with  respect 
to  the  direction  of  flow  by  a  mark  that  had  previously  been  made  on  the  end  of  the  rod. 

The  test  was  run  for  a  sufficient  length  of  time,  from  3/4  to  2  hr,  to  remove  from  0.01  in. 
to  0.02  in.  from  the  surface  of  the  rod.  Temperature,  pressure,  and  flow-rate  measurements 
were  taken  several  times  during  this  Interval.  At  the  end  of  the  test,  the  rod  was  removed 
and  the  end  tape  stripped  off.  A  paper  disk,  marked  in  10-increments,  was  glued  onto  one 
end  with  the  zero  mark  pointing  toward  the  oncoming  air  stream.  The  depth  of  surface  re¬ 
moval  was  measured,  with  reference  to  the  unexposed  surface  at  the  ends,  at  10-deg  Intervals. 
These  measurements  were  made  by  use  of  a  comparator  which  projected  an  enlarged  view  of  the 
rod  on  a  screen.  These  measurements  were  made  at  the  center  cross  section  but  were  checked 
by  measurements  near  the  end.  There  was  no  significant  variation  between  measurements  at 
the  two  points. 

Evaluation  of  Data 


The  local  mass -transfer  coefficients  were  calculated  from  the  equation 


y  /^n 

n»  M 
*  n 


(8) 
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The  average  coefficients  were  calculated  in  the  same  way  using  the  average  depth  of 
naphthalene  removal  as  computed  by  a  planimeter  Integration  of  the  local  values.  The 
corresponding  J -factors  were  then  found  from  the  defining  equation,  using  experimental 
values  for  flow  rate,  and  air  and  naphthalene  properties  estimated  from  tempera ture 
and  pressure  data.  For  comparison  with  tests  performed  by  other  workers,  the  dimen¬ 
sionless  variable  Nu/Pr  was  calculated  from  the  relation 

Nu/Pr  ^  «*  J.  Re 

n 


taking  Jh  =  Jm.  Heat-transfer  coefficients  also  may  be  calculated  from  Jh  by  means 
of  the  defining  equation.  The  values  for  the  molecular  diffusion  coefficient  were  com¬ 
puted  from  the  formula  given  by  Sherwood  and  Pigford  (jj) .  The  vapor  pressures  were 
calculated  from  the  correlation  given  in  "International  Critical  Tables"  (J) • 

Discussion  of  Results 


Average  heat-transfer  coefficients  as  measured  by  the  mass -transfer  technique 
are  plotted  in  dimensionless  form  in  Fig.  6  and  compared  with  those  obtained  by  transient 
tests.  It  Is  seen  that  the  mass-transfer  results  are  subject  to  more  scatter  than  the 
transient  results  and  average  about  20  per  cent  higher.  Comparison  is  also  made  with 
the  Colburn  equation  and  with  the  results  of  Winding  and  Cheney  (8) .  The  Winding  and 
Cheney  experiments  also  made  use  of  the  naphthalene  mass-transfer  technique  but  were  per¬ 
formed  on  a  bank  of  tubes  1.5  in.  diam  whose  center-to-center  spacing  within  a  row  was 
2.25  in.,  with  rows  3.00  in.  apart.  Winding  and  Cheney  reported  the  correlation 

Nu  =  0.55  (Re)0-55 


which  becomes 


0.55  1/3 

Nu  =  0.622  (Re)  (Pr) 


assuming  Pr  »  0.70.  The  results  reported  here  lie  also  above  both  of  these  correlations. 
There  are  several  factors  that  the  authors  believe  tend  to  make  it  more  difficult  to  obtain 
accurate  information  by  the  mass-transfer  technique.  These  are  as  follows* 

1  A  slight  darkening  of  the  naphthalene  during  the  casting  process  indicates  the 
possibility  of  chemical  change  that  might  Influence  the  physical  properties. 

2  Although  the  rods  were  smooth  immediately  after  casting,  varying  degrees  of 
surface  roughness  developed  during  sublimation. 

3  The  vapor  pressure  of  naphthalene  changes  rapidly  with  temperature.  Data  on 
vapor  pressure  in  the  range  of  room  temperatures,  and  below,  are  meager. 

4  Some  of  the  rods  tended  to  be  porous  or  flaky  after  casting.  Such  rods,  in 
general,  were  rejected,  but  detection  of  flaws  was  not  always  possible. 

In  spite  of  the  errors  in  absolute  values  that  appeared  to  creep  Into  these  ex¬ 
periments,  the  mass-transfer  technique  is  nevertheless  believed  by  the  authors  to  be  one  of 
the  better  methods  of  detecting  normalized  variations  in  heat-transfer  coefficients.  Plots 


of  local  heat-transfer  coefficients  for  row  3  are  shown  In  Pigs.  7  and  8.  Similar  plots 
for  rows  5  and  7  are  shown  In  Figs.  9  and  10.  The  maximum  deviations  from  the  average 
coefficients  are  given  In  Table  1. 


Table  1  Maximum  Deviations  from  Average  Coefficients 
Row _  _ Re _  Max,  Deviation  from  Ave.  h  (Per  Cent) 


3 

36000 

-  43800 

+79-5 

-52.2 

3 

73700 

+67.0 

-47.3 

5  and  7 

36000 

-  43800 

+64.4 

-54.5 

5  and  7 

73700 

+57.8 

-39-2 

Point-to-point  variations  are  somewhat  greater  than  those  reported  by  Winding  and 
Cheney,  but  are  in  qualitative  agreement.  The  results  are  not,  however.  In  agreement  with 
the  findings  of  Dwyer,  Sheehan,  and  Welsman  .  Dwyer,  et  al.  measured  local  coefficients 
from  staggered  tubes  to  water  at  a  temperature  of  360  P  using  the  same  basic  equipment  as 
is  described  in  (2) .  Maximum  deviations  from  the  average  aoefflalent 
were  reported  to  be  +9  per  cent  and  -6.1  per  cent  Re  *  27,600,  and  +6  per  cent  and  -1.4  per 
aent  at  Re  ■=  74,200.  At  the  same  time,  the  sharp  variation  that  the  authors  found,  between 

80  and  140  P  from  the  forward  stagnation  point,  were  not  reported  by  Dwyer,  et  al. 
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Fig.  2  Drawing  of  thermal  capacitor  for  transient  tests 
Fig.  3  Temper  attire -time  history  of  test  rod 

Fig.  4  Transient  test  results  for  third  row 

Fig.  5  Transient  test  results  for  fifth  and  seventh  rows 

Fig.  6  Average  heat  transfer  measured  by  mass-transfer  method 
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